A radio-frequency sheath model for complex waveforms 
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Abstract 

Plasma sheaths driven by radio-frequency voltages occur frequently, in contexts ranging from 
plasma processing applications to magnetically confined fusion experiments. These sheaths are cru- 
cial because they dominantly affect impedance, power absorption, ion acceleration and sometimes 
the stability of the nearby plasma. An analytical understanding of sheath behavior is therefore 
important, both intrinsically and as an element in more elaborate theoretical structures. In prac- 
tice, these radio-frequency sheaths are commonly excited by highly anharmonic waveforms, but 
no analytical model exists for this general case. In this letter we present a mathematically simple 
sheath model that can be solved for essentially arbitrary excitation waveforms. We show that this 
model is in good agreement with earlier models for single frequency excitation, and we show by 
example how to develop a solution for a complex wave form. This solution is in good agreement 
with simulation data. This simple and accurate model is likely to have wide application. 
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In many radio-frequency discharges the sheath is the most important region, because 
impedance, power absorption and ion acceleration are dominated by sheath processes [l, 2|. 
There are other contexts where radio-frequency sheath physics is a concern, for example 
when understanding the physics of heating in certain fusion plasmas [3 5 1 - Consequently, 
models of the sheath are important, either in themselves or as elements in more complex 
situations, over a broad area of plasma physics. The problem has been considered on a 
number of occasions, e.g. jcj-jj. Analytical models are particularly useful for developing 
physical insight and expressing the relationships between parameters in a clear way, but 



such models have proved elusive. Liebermanjl, 



9[ supplied an analytical model for a radio- 



frequency sheath driven by a single frequency, but in practice much more complex waveforms 



frequently occur 



1CH14| . There has been limited success in generalizing the Lieberman model 



15Ml7j]. So there is essentially 



to cover these cases, because of mathematical complexities 
no sheath model available to describe many modern experiments. In this paper we present 
a new analytical sheath model, based on a simpler mathematical framework than that of 



Lieberman 



LyticaJ 
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For the single frequency case, this model yields scaling laws that are 
identical in form to those of Liebermanjl, 9], differing only by numerical coefficients close to 
one. However, the new model may be straightforwardly solved for almost arbitrary current 
waveforms, and may be used to derive scaling laws for such cases. 

Fig. [1] is a schematic representation of the charged particle densities and fields that occur 
in a radio-frequency sheath. For a sheath in this regime, the ion motion is determined by 
the time averaged field, while the electrons respond to the instantaneous field. A model 
describing such a sheath therefore has a time-averaged part and a time- dependent part, 
which must be consistent. We therefore insist on the same maximum sheath width, s m , 
in both cases. The principal parameter determining s m in the time-averaged sense is the 
time-averaged sheath voltage, V, while in the time-dependent model, s m is a function of V , 
the maximum sheath voltage. The simplest way to satisfy the constraint is to choose 



T T S 1 

V po 



(1) 



where p is the time-averaged charge density, p is the charge density when the sheath voltage 
is Vo, and £ emerges as a key parameter of our model. We note that p 7^ P because of the 
time dependence of the electron density in the sheath region. Since the ion density rij is 
time- independent, these relations imply n e = (1 — which is our central approximation. 
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We now consider specifically a sheath adjacent to a plasma of density uq and electron 
temperature To. Ions flow into the sheath at x = and are absorbed at an electrode at 
x = s m . The constant ion current density is J, = en UB, where m b = ^Jk^T^jM and M is 
the ion mass. Now the governing equations for the time-averaged ion motion are 



n u B 
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dx 2 e e 

where Ui is the (time independent) ion drift velocity and is the time-averaged potential, 
and we assume that eV ^> k B T . These equations are those of the Child-Langmuir sheath 
model [l[ [2I, [isl isl, with the addition of the parameter £, and we can at once write down 
the solutions 
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where the boundary condition E(x = 0) = has been used, and where Ki = 4/(9£). The 
time- dependent field and potential can now be determined by integrating Poisson's equation 
again with the assumption that 

r if s < x < s m ; 

n e = \ 

{ Hi otherwise, 

where s(t) is the position of the sheath edge, at which point E = and ( 
for s < x < s m '- 



(9) 

0. We obtain 
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so that the time dependent sheath voltage is 

4 / s 



V(t) = V 
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Now 



so that 
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Since < s/s m < 1, once J(i) is chosen, s(t)/s m is fully defined and we can express 
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(14) 
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Hence, once J(t) has been given, all the remaining quantities can be calculated without 
further assumption or approximation, as the following example shows. 

For the single frequency case treated by Lieberman [9j we choose J(t) = — Jo sin cut. From 
Eqs. (USD and (TUD we find 



s(t) 
Jo 



1 — cos cot) 3 
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(16) 
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where K cap = 4/3. The voltage waveform is obtained by inserting Eq. (fl6l) into Eq. (Tl2 
Combining Eqs. (jSJ) and ( TT7|) gives the sheath maximum expansion as a function of Jo, 



i^J 3 



(19) 



ee fcijToCU 3 ng' 

where -fT, = 4^/3. These expressions are identical with those of Lieberman , apart 
from numerical coefficients close to unity. Table [T] compares £, Ki, K ca _ p and K s for the two 
models and shows that they are not significantly different. Similarly, the voltage waveforms 
are almost identical in the two models. This may seem surprising, in view of the apparently 
bold approximation of Eq. ([T]). However, Brinkmann 20] has shown that the approximation 
of Eq. (jnD is important, and because this approximation is integral to the Lieberman model, 
the accuracy of that model is not significantly better than the present one. However, the 
present model can be solved for a far greater range of waveforms. In particular, the current 
density can be expressed as an arbitrary Fourier series, leading, for example, to models for 
multiple-frequency excitation that are free of inconvenient restrictions on the component 
amplitudes (such as occur in dual-frequency generalizations of Lieberman's model 15Hl7j]). 
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For example, we can choose 



J(t) 



-Jo sinwot — J\ smuiit 



(20) 



and find at once 
a(t) 
Jo 
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where we note that eq. [24] yields a result never different by more than 10 % from the 
single frequency result, which value can therefore be used for all practical purposes. These 
formulae are therefore more generally applicable, less cumbersome and obtained with less 



mathematical exertion than those previously given 



15ril7| 



As a third example we consider a sheath excited by the pulsed waveform 

J(t) 
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10Nl2|. We assume that this pulse 



which is representative of several topical experiments 
is repeated at intervals t p <C t w , such that successive pulses do not appreciably overlap. In 
this case we find 
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The Child law, and therefore the sheath width, are found by inserting Eq. (|28|) into Eq. ([5]). 
We have investigated the utility of this model by comparison with particle-in-cell simulation 
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data 



21 



22J. These simulations treated a plasma formed in a space between two plane 



parallel electrodes separated by 6.7 cm, filled with argon gas at a pressure of 10 mTorr, and 
excited by a current density with the form of Eq. (I2"6"j) . The peak current density ranged 
from 5 to 70 A m -2 and the pulse width from approximately 1 to 10 ns. These conditions 
lead to Tq « 1.5 eV and Uq ~ 3 x 10 14 -3 x 10 15 m -3 . In figs. [2J [3] and H] we compare these 
simulation results with the predictions of the present model, and we find good agreement, 
both for time time dependent currents and voltages, and scaling laws. 

To summarize, in this paper we have developed a new sheath model that can be ex- 
pressed in a small number of straightforward equations. For the single-frequency case, the 
model agrees well with the Lieberman model. The Lieberman model, in spite of its ele- 
gant construction, is mathematically complex and has proved resistant to generalization. 
The present model, however, is readily adaptable to a wide range of complex excitation 
waveforms, and leads to results in good agreement with simulations. Thus this model pro- 
vides a valuable tool for designing and understanding experiments involving non-sinusoidal 
excitation of radio-frequency sheaths. 
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FIG. 1: Upper panel: Charged particle densities, showing the ion density, rii (solid line) the time 
averaged electron density, n e (dotted line) , and the electron density, n e at the instant when s/s m = 
0.3 (dashed line). Lower panel: Electrostatic potential, showing the time averaged potential, (f> 
(solid line) and the instantaneous potential when s/s m = 0.3 (dashed line). 
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FIG. 2: Comparison of simulation results (solid lines) with the analytical theory of the text (dashed 
lines) for sheath voltage (upper panel) and sheath current density (lower panel). For this case 
t w = 5.2 ns. The electrical control parameter for the sheath model is V, the time averaged sheath 
voltage, which is here chosen to be the same as in the simulation. 
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FIG. 3: Comparison of the maximum radio- frequency current density found in simulation (hor- 
izontal axis) with the result computed from eq. [29] (vertical axis). The solid line denotes ideal 
agreement between theory and simulation. 
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FIG. 4: Comparison of the ion current density found in simulation (horizontal axis) with the result 
computed from eq. [5] and eq. [28] (vertical axis). The solid line denotes ideal agreement between 
theory and simulation. 
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